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Abstract 

Antibody-based technology is the main method for diagnosis and treatment of snake bite envenoming currently. However, 
the development of an antibody, polyclonal or monoclonal, is a complicated and costly procedure. Aptamers are single 
stranded oligonucleotides that recognize specific targets such as proteins and have shown great potential over the years as 
diagnostic and therapeutic agents. In contrast to antibodies, aptamers can be selected in vitro without immunization of 
animals, and synthesized chemically with extreme accuracy, low cost and high degree of purity. In this study we firstly 
report on the identification of DNA aptamers that bind to p-bungarotoxin ((3-BuTx), a neurotoxin from the venom of 
Bungarus multicinctus. A plate-SELEX method was used for the selection of P-BuTx specific aptamers. After 10 rounds of 
selection, four aptamer candidates were obtained, with the dissociation constant ranged from 65.9 nM to 995 nM measured 
by fluorescence spectroscopy. Competitive binding assays using both the fluorescently labeled and unlabeled aptamers 
revealed that the four aptamers bound to the same binding site of P-BuTx. The best binder, PB-1, bound specifically to p- 
BuTx, but not to BSA, casein or a-Bungarotoxin. Moreover, electrophoretic mobility shift assay and enzyme-linked aptamer 
assay demonstrated that PB-1 could discriminate B. multicinctus venom from other snake venoms tested. The results 
suggest that aptamer pB-1 can serve as a useful tool for the design and development of drugs and diagnostic tests for p- 
BuTx poisoning and B. multicinctus bites. 
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Introduction 

Snakebite is a major public health issue that affects residents of 
rural communities living in tropical and subtropical countries [1]. 
It was estimated that, globally, between 1 .2 million and 5.5 million 
people are bitten by snakes annually, resulting in about 20,000- 
125,000 deaths, and 400,000 permanent disabilities [2-4]. The 
administration of animal-derived antivenoms is the mainstay 
treatment of snakebite envenoming. Conventional antivenoms are 
prepared from sera of large animals, usually horses or sheep, 
hyperimmunized with relevant snake venoms. After collection of 
blood or plasma, the plasma is fractionated to extract and purify 
the active immunoglobulin or immunoglobulin fragments (F(ab')2), 
and then the antivenom is freeze-dried for conservation and 
transportation. Due to its complex and time consuming process of 
production, antivenom is so expensive that many people actually 
cannot afford the high costs. Moreover, antivenom therapy of 
snakebite has long been associated with a high incidence of 
adverse reactions such as pruritus, urticaria and potentially fatal 
anaphylaxis [5,6]. The complexity of the production of anti- 
venoms, the decreasing number of producers and the fragility of 
the production systems in developing countries jeopardize the 



availability of effective antivenoms in Africa, Asia, the Middle East 
and South America [7]. 

Identification of snakebite is also difficult as lack of a rapid and 
reliable diagnostic test. Polyvalent antivenom is often used when 
the snakebites can not be identified, theoretically increasing the 
risk of late antivenom reactions because a higher dose of foreign 
protein is administered [5]. Current efforts in the diagnosis of 
snake envenomation have also been mainly focused on the 
development of antibody-based immunoassays. Snake venoms are 
highly complex mixtures that tend to have many common antigens 
among snakes from the same or even different genera [8,9]. Direct 
use of crude antibody in immunoassays may have low specificity in 
species detection of snakebites. Monoclonal antibodies against a 
single venom component or species-specific polyclonal antibodies 
produced by passing the crude antisera through the medium 
coupled with heterologous venoms were used to reduce cross 
reactivities [8,10,11]. 

A new class of short single-stranded oligonucleotides (RNA or 
DNA) termed aptamers has quickly emerged as a novel and 
powerful class of ligands with excellent potential for diagnostic and 
therapeutic applications [12,13]. Aptamers are selected in vitro 
through the systematic evolution of ligands by exponential 
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enrichment (SELEX) process [14,15], that have high affinity and 
specificity to a wide range of targets, such as small molecules, 
proteins, or whole cells [16-18], As termed "chemical antibod- 
ies", aptamers have several advantages over their counterparts. 
They can be selected in vitro from random DNA or RNA 
libraries without immunization of animals, and synthesized 
chemically in a readily scalable process with extreme accuracy, 
low cost and high degree of purity. They display low to no 
immunogenicity even administered in high doses. Aptamers can 
also be easily modified by various dyes or functional groups 
during chemical synthesis and immobilized on transducing device 
widely used in biosensors. 

The many-branded krait Bungarus multicinctus is widely 
distributed throughout Southeast Asia. The snake possesses 
extremely toxic venom, one of the most potent of any land snake. 
Bitten by B. multicinctus does not give rise to swelling or necrosis 
at the site of the bite, but cause severe neuromuscular blockade, 
resulting in respiratory failure and fatality [19-21]. B. multicinctus 
was responsible for 8.12% of snakebites and the highest mortality 
in China [22]. The venom of B. multicinctus consists of both pre- 
and postsynaptic neurotoxins, including a-, ($-, y- and K- 
bungarotoxin. Although oe-bungarotoxin (a-BuTx) is the major 
component (61%) of the B. multicinctus venom [23], fi-bungaro- 
toxin (fi-BuTx) is more important and lethal than a-BuTx from a 
toxicological or clinical perspective [24]. (3-BuTx constitutes > 
20% of the protein content of the crude venom [25]. The 
minimum lethal dose of P-BuTx by intraperitoneal injection of 
mice is 0.01 ug/g [26], while that of a-BuTx is 0.139 ug/g [27]. 
(3-BuTx is a presynaptic neurotoxin with a molecular weight of 
21.8 kDa, consisting of two dissimilar polypeptide chains, the 
phospholipase A2 subunit A chain (~ 1 4 kDa) and the non- 
phospholipase A2 subunit B chain (~7 kDa), cross-linked by an 
interchain disulfide bond [26,28,29]. |3-BuTx presynaptically 
inhibits acetylcholine release in both the peripheral and central 
nervous systems [30]. 

Until now, several aptamers targeting biotoxins have been 
selected [31], including bacterial toxins [32], mycotoxins [33] and 
plant toxins [34]. However, few aptamers are developed against 
snake venom components. The only aptamer against snake venom 
component developed to the present is a DNA aptamer against ot- 
BuTx [35]. But the binding affinity of this aptamer is relatively 
low, with a dissociation constant of 7.58 U.M, thus may limit its use 
in biological applications. Here in the present study we aim at 
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SELEX round 

Figure 1. Binding affinity of ssDNA pools of each round of 
selection. Ten nM of ssDNA pool was incubated with 100 nM of (S- 
BuTx and the ratio of ssDNA bound to p-BuTx was determined by 
nitrocellulose filter binding assay. 
doi:1 0.1 371 /journal.pone.01 05404.g001 



generating aptamers with high affinity and specificity against (3- 
BuTx that can be applied to the diagnosis and/or therapy of 
snakebites. By using a plate-SELEX process, aptamers were 
selected from a random ssDNA library. The dissociation constant 
values were determined by fluorescence anisotropy measurements 
and found to be in the range of 10 M. The aptamers bind with 
high affinity only to fi-BuTx and fi-BuTx containing B. multi- 
cinctus venom, but not to a-BuTx or other snake venoms tested. 
The ability of the aptamers to specifically recognize B. multi- 
cinctus venom offering the potential for application to the 
diagnosis of B. multicinctus bites. 

Materials and Methods 

Reagengts 

The 80-nt single-stranded DNA (ssDNA) library contained a 
randomized region of 40-nt central region flanked by two 
conserved 20-nt primer hybridization regions in each end (5'- 
AGCAGCACAGAGGTCAGATG-N40-CCTATGCGTGCTA- 
CCGTGAA-3'). The fluorescein labeled forward primer (F-fam: 
5 '-FAM-AGCAGCACAGAGGTCAGATG-3 ') and the phos- 
phate labeled reverse primer (R-pho: 5'-phosphate-TTCACGG- 
TAGCACGCATAGG-3') were used in PCR to obtain double- 
stranded DNA (dsDNA). The library and all primers were 
synthesized by Sangon Biotech (Shanghai, China). Lyophilized 
crude venoms were obtained from a serpentarium from Jiangxi 
Province, China. BSA and casein were purchased from the 
Sigma-Aldrich Co. (USA). (5-BuTx was from Santa Cruz (USA), 
and a-BuTx was from Abeam (UK). Premix Taq polymerase, 
pMD19-T vector, and DH5a competent cells were purchased 
from TaKaRa Biotechnology Co. Ltd. (Dalian, China). All the 
other reagents used for chemical and biological characterization 
were of analytical grade. 

SELEX Procedure 

To select aptamers directed to [5-BuTx, 10 rounds of plate 
SELEX were carried out. A 96-well ELISA polystyrene plate 
(Corning, USA) was coated at 4°C overnight with 100 uL of 
10 ug/mL (3-BuTx in 0.05 M carbonate-bicarbonate buffer 
(1.59 g/L Na 2 C0 3 and 2.93 g/L NaHCO :j , pH 9.6). After 
coated, the wells were blocked at room temperature (RT) for 
2 h using binding buffer (20 mM HEPES, pH 7.4, 150 mM NaCl, 
5 mM KC1, 2 mM MgCl 2 and 2 mM CaCl 2 ) containing 1 mg/ 
mL BSA for rounds 1-5. For rounds 6-10, wells were blocked with 
binding buffer containing 1 mg/mL casein. This switch of 
blocking agent was intended to decrease background binding of 
aptamer pools to the plate surface [36]. The ssDNA library or 
aptamer pool in binding buffer was heated at 90°C for 10 min in a 
thermal cycler, immediately left to cool on ice for another 1 0 min, 
and then kept at RT for 10 min. For the first round of selection, 
1 nmol of ssDNA library (~10 14 different sequences) was 
incubated with |3-BuTx coated wells at RT for 1 h. For the next 
9 rounds of selection, ssDNA pool used for incubation was 
gradually decreased to 50 pmol. After incubation, the supernatant 
was discarded and the wells were washed 6 times with 200 uL of 
washing buffer (binding buffer containing 0.05% Tween 20), and 
further washed 3 times with 200 u.L of ddH 2 0. To elute the 
bound ssDNA, 1 00 u.L of ddH 2 0 was added to the wells (each 
well was separated from the strip and sealed with silver paper) and 
placed in a Thermomixer Compact (Eppendorf, Germany) at 
temperature of 95°C for 10 min. The supernatant was collected 
and served as the template for the amplification of fS-BuTx-bound 
ssDNA. 
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Table 1. Sequences of selected aptamers and their dissociation constants (K d ). 




Name Sequence of variable region 


Percentage of total 
sequences(%) 


K d InM) 


PB-1 lGTTTtTCCCdTTGTC|GCTrnlGGTllCGTTCTGCCTCTATCT 40 nt 


73.0 


65.9 


PB-20 ATAOTC ATiil 1 1 CiTI I ( i l( TGGC 1 1 1 l|l GiGGTTTGTGCAGTATTATGAAC 51 nt 


8.1 


83.8 


PB-1 9 TTTGGTGTGGATCCTGAACATTTATATTCTTTCG 1 1 1 1 1 1 40 nt 


16.2 


530 


PB-32 GCAATGCACCTTTGTCTCTTATAGTTTATT TTTTGCCTT 39 nt 


2.7 


995 



The boxed sequences were constant regions between aptamers (3B-1 and PB-20. 



During every round of selection, a pilot PCR was first 
performed to determine the optimal cycle number of amplification 
as described by Lou et al. [37] . Then dsDNA was amplified under 
the optimal cycles. One hundred microliters of the PCR mixture 
contained 50 uL Premix Taq solution, 4 |jL template, 1 uL of 
10 uM fluorescein labeled forward primer F-fam, 1 uL of 10 uM 
phosphate labeled reverse primer R-pho, and 44 ilL of ddH 2 0. 
PCRs conditions were as follows: preincubation at 94°C for 5 min; 
optical number of cycles consisting of denaturation at 94°C for 
0.5 min, annealing at 56°C for 0.5 min and extension at 72°C for 
0.5 min; and final extension at 72°C for 10 min. The dsDNA 
PCR product was digested by X exonuclease (NEB, UK) for 2 h at 
37°C to generate ssDNA, according to the manufacturer's 
protocol. Reaction was stopped by 10 min incubation at 75°C. 
ssDNA was precipitated with dehydrated alcohol and washed with 
70% alcohol. After centrifugation, the pellet was dissolved in 
binding buffer and used as the enriched library for the next 
selection round. The concentration of ssDNA was evaluated by 
electrophoretic analysis [38]. Each sample (10 uL) was compared 
with 10 of ssDNA library of known concentration (20, 40, 60, 
80, 100 nM) using agarose gel electrophoresis. Electrophoretic 
analysis of the band intensities was performed using Image J 
program. 

After 10 rounds of selection, the collected aptamer pool was 
PCR amplified with unlabeled forward and reverse primers. The 
dsDNA product was subjected to 3% agarose gel and the 
corresponding strand was excised and purified using the QIAEX 



II Gel Extraction kit (Qiagen, Germany) according to the 
manufacturer's instructions. Purified PCR product was ligated 
into the pMD19-T vector and transformed to DH5ot competent 
cells by heat shock. Forty-five colonies were randomly picked and 
sequenced by Sangon Biotech (Shanghai, China). Secondary 
structure analysis was performed with Mfold software (http:// 
mfold.rit.albany.edu/?q = mfold/DNA-Folding-Form). 

Nitrocellulose filter binding assay 

The 0.22-u.m HAWP nitrocellulose filters (Millipore) were pre— 
treated with alkali to reduce the non-specific adsorption of nucleic 
acids as described [39] . The filters were soaked in 0.5 M KOH at 
RT for 20 min, then washed extensively with ddFLO. After 
washing for 45 min in binding buffer by shaking, the filters were 
stored in new buffer at 4°C. 

For measuring the aptamer pool binding, 1 mL of 10 nM 
fluorescein-labelled aptamer pool in binding buffer was heated to 
90°C for 10 min and rapidly put on ice for 10 min, then incubated 
10 min at RT. P-BuTx was then added to a final concentration of 
100 nM and incubated at RT for 40 min. Controls did not 
contain P-BuTx. Afterwards, the mixture was vacuum filtered over 
a pretreated nitrocellulose, and washed with 1 mL of binding 
buffer. Total volume of 2 mL of aptamer DNA eluted from the 
filter was determined by fluorescence measurement using a Perkin- 
Elmer LS55 Fluorescence Spectrophotometer. The excitation 
wavelength was set at 492 nm and the emission wavelength at 
5 1 8 nm. Slits for both the excitation and the emission were set at 
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Figure 2. Dissociation constant (K d ) measurements by fluorescence anisotropy. A) (-SB-1; B) pB-20. Two nM of fluorescein labeled aptamers 
were titrated with increasing concentrations of p-BuTx. The anisotropy change (delta r) was the anisotropy of the initial fluorescein labeled aptamer 
subtracted from the anisotropy value at each P-BuTx concentration. Non4inear regression analysis revealed the K d of pB-1 was 65.9±7.8 nM, and the 
K d of PB-20 was 83.8±4.3 nM. Error bars represent standard deviations from three repeated experiments. 
doi:1 0.1 371 /journal.pone.01 05404.g002 
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Figure 3. Determination of the binding sites of selected 
aptamers by competitive anisotropy assay. Two nM of fluores- 
cently labeled aptamer (|3B-1, (3B-20, pB-19 or (3B-32) was firstly 
incubated with 288 nM of (3-BuTx, and then titrated with increasing 
concentrations of unlabeled aptamer (3B-1. Labeled PB-1 was also 
titrated with unlabeled random ssDNA as control. 
doi:1 0.1 371 /journal.pone.01 05404.g003 

10 nm, and the integration time was set at 10 sec. The 
fluorescence of sample elution subtracted from control elution 
gave the fluorescence of P-BuTx-bound ssDNA that maintained 
on the filter. 

To measure the aptamer binding to other relative proteins, 
fluorescein labeled aptamer was chemically synthesized after 
sequencing. Ten nM of fluorescein labeled aptamer was mixed 
with 100 nM of P-BuTx, BSA, a-BuTx and 2.1 ug/mL of casein 
in binding buffer, respectively, and then subjected to nitrocellulose 
filter binding assay using the same above procedure. 

Fluorescence anisotropy measurements 

Fluorescence anisotropy measurements were performed on 
Perkin-Elmer LS55 Fluorescence Spectrophotometer with excita- 
tion and emission polarizers in L-format, and the settings were the 
same as above. The G factor was automatically calculated by the 
instrument. All experiments were carried out at RT. The 
concentration of fluorescein labeled aptamers was 2 nM (in 



binding buffer) and the total volume was 2 mL. Different 
concentrations of (5-BuTx were titrated into the aptamer solution 
and incubated for 2 min, and then the fluorescence anisotropy 
values were measured. Four anisotropy measurements were taken 
each time, and the experiments were repeated three times. The 
relative standard deviation was <3% for all the measurements. 
The change in anisotropy was the average anisotropy of the initial 
fluorescein labeled aptamer subtracted from the average anisot- 
ropy value at each (5-BuTx concentration. The dissociation 
constant (K<j) was determined by non-linear regression analysis 
with software SigmaPlot 10.0 by using a single-site ligand binding 
model. 

For competitive displacement experiments, 2 nM of fluorescein 
labeled aptamer (PB-1, PB-20, PB-19 or PB-32) and 288 nM of P- 
BuTx were premixed in 2 mL of binding buffer and following that, 
titration experiments were conducted with increasing amount of 
unlabeled PB-1 respectively. Fluorescence anisotropy values were 
measured as mentioned above. As control, the curve generated by 
fluorescein-labeled PB-1 against unlabeled starting ssDNA library 
was also measured. 

For specificity assays, the anisotropy of 2 nM of fluorescein 
labeled aptamer PB- 1 in 2 mL of binding buffer was measured. 
And then, 50 nM of P-BuTx or 500 nM of BSA and a-BuTx, or 
1 0 |Xg/mL of casein was added, respectively. The corresponding 
changes of anisotropy were recorded. 

Electrophoretic Mobility Shift Assay (EMSA) 

Crude venoms of B. multicinctus , B. fasciatus, Deinagkistrodon 
acutus, Gloydius brevicaudus and Naja Naja Atra were resolved in 
binding buffer, and centrifuged at 12 000 g for 10 min, and the 
supernatant was collected and stored at — 20°C. Fluorescein 
labeled aptamer (3B-1 (100 nM) was mixed with P-BuTx (240, 480 
or 960 nM) or B. multicinctus venom (15, 30 or 50 (Xg/mL) in 
binding buffer of 20 uL total volume, respectively. The binding 
reaction was carried out at RT for 40 min and then analyzed by 
1% agarose gel containing 0.5 trig/ml of ethidium bromide. 
Binding reaction was also performed by using 100 nM of PB-1 and 
50 ng/mL of the five snake venoms mentioned above. 

Enzyme-linked aptamer assay (ELAA) 

Aptamer DNA was labeled with biotin during chemical 
synthesis. BSA, casein, a-BuTx and different crude snake venoms 
were diluted to 10 u,g/mL in 0.05 M carbonate-bicarbonate 



B-BuTx 



fJ-BuTx: 100 nM (2.1 jig/mL) 
a-BuTx: 100 nM 
BSA: 100 nM 
Casein: 2.1 fig/mL 



B 0.025 



P-BuTx: 50 nM (1.05 ug/mL) 
a-BuTx: 500 nM 
BSA: 500 nM 
Casein: 10 ug/mL 



Figure 4. Binding affinity of the aptamer pB-1 to related proteins. A) Ten nM of fluorescein labeled (3B-1 was mixed with 1 00 nM of P-BuTx, 
BSA, oe-BuTx and 2.1 ug/ml_ of casein, respectively, and then nitrocellulose filter binding assay was performed to calculate the amount of pB-1 bound 
to tested protein; B) Two nM of fluorescently labeled aptamer (3B-1 was mixed with 50 nM of (3-BuTx and 500 nM of oi-BuTx and BSA, and 10 |rg/ml_ of 
casein, respectively, the anisotropy changes of (3B-1 were determined. Values represent means ± standard errors. 
doi:1 0.1 371 /journal.pone.01 05404.g004 
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Figure 5. EMSA analysis of the binding of aptamer pB-1 to p-BuTx and five snake venoms. One hundred nM of |>B-1 was incubated at RT 
for 40 min with 240 nM (2), 480 nM (3), 960 nM (4) of P-BuTx, 1 5 ng/ml_ (5), 30 ug/mL (6), 50 ug/mL (7 and 9) of fi. multicinctus venom, and 50 ug/mL 
of venom of 6. fasciatus (10), Naja Naja Atra (1 1), D. acutus (12), and G. brevicaudus (13). Lane 1 and Lane 8: no protein or venom added. F: free DNA; C: 
DNA/protein complex; S: DNA stuck in well. 
doi:1 0.1 371 /journal.pone.01 05404.g005 



buffer, pH 9.6, and incubated in an ELISA polystyrene plate strip 
(Corning, USA) overnight at 4°C to allow plate binding. Plate 
wells were washed three times with 200 uL of washing buffer, and 
then blocked at RT for 2 h with 100 uL of 1 mg/mL BSA in 
binding buffer. After washing three times in washing buffer, 
1 00 U.L of 1 0 pM biotin-labeled aptamer was added to individual 
wells and incubated at RT for 40 min. After the above procedures, 
individual wells were washed three times with washing buffer to 
remove unbound ssDNA. Then 1 00 |J.L of a 1:1 000 dilution of 
horseradish peroxidase (HRP)-labeled Streptavidin (Solarbio, 
China) in binding buffer was added to the wells. Following a 40- 
min incubation at RT, the wells were washed three times and 
1 00 uL of freshly prepared substrate OPD solution (Sigma, USA) 
was added according to the manufacturer's instruction. The wells 
were kept in dark at RT for 15 min and 50 uL of 2 M H2SO4 was 
added to each well to stop the reaction. 

Results and Discussion 

SELEX Selection 

Starting from a library of ~10 14 different sequences of 
oligodeoxyribonucleotides, we selected aptamers binding to (3- 
BuTx, using plate SELEX technique. The critical step in SELEX 
is the partitioning of the target-bound oligonucleotides from the 
unbound ones. There are several ways to carry out the separation 
process, such as using of nitrocellulose membrane [40], affinity 
column [41], magnetic beads [42] or capillary electrophoresis 
[43]. Immobilization of target protein on a 96-well microtiter plate 



provides a simple and rapid way to separate aptamer-target 
complex from free DNAs. We first adsorbed (3-BuTx to 
polystyrene 96-well plates, and then incubated with the ssDNA 
library pool containing a 40-nucleotide random sequence region. 
Following washes by washing buffer and ddH 2 0 to remove 
unbound ssDNA, the P-BuTx-bound aptamers were eluted by 
heat shock at 95°C with ddH z O. High concentration of 
guanidinium isothiocyanate buffer was often used in elution of 
aptamers from the plate, and a phenol-chloroform extraction and 
an alcohol precipitation procedure was performed to purify the 
ssDNA [44,45]. We found that a heat shock with ddH 2 0 was 
sufficient to elute bound ssDNA from the plate, and it could be 
directly served as the template in the following PCR reactions 
without any purification steps. There was probably some residual 
(3-BuTx in the elution solution, but the [3-BuTx was denatured, so 
we thought it would have no significant interference in the PCR 
amplifications. 

The selected ssDNA pool was amplified using pilot PCR in 
ranges of 14 to 26 PCR cycles in order to optimize the 
amplification cycle number under which neither shorter nor 
longer PCR byproducts appeared (Figure SI). The forward and 
reverse primers were fluorescein and phosphate labeled at 5' end, 
respectively. Then the dsDNA PCR product was amplified under 
the optimal cycles and subjected to X exonuclease digestion to 
produce fluorescein labeled ssDNA. Phosphorylated strands of 
DNA are selectively digested by the enzyme, which has gready 
reduced activity on ssDNA and non-phosphorylated DNA. The 
generated ssDNA was used for the next cycle of SELEX. The 



Figure 6. Identification of p-BuTx and B. multicinctus venom using aptamer pB-1 via ELAA. Snake venoms and proteins of 10 ug/mL were 
coated onto ELISA plate wells and incubated with 10 pM biotin-labeled aptamer pB-1, respectively. Then HRP-labeled Streptavidin was added and 
substrate OPD was used to achieve colorimetric responses. Control well was coated with p-BuTx but no aptamer was added. Well 1 : P-BuTx; Well 2: ot- 
BuTx; Well 3: BSA; Well 4: casein; Well 5: B. multicinctus; Well 6: B. fasciatus; Well 7: D. acutus; Well 8: G. brevicaudus; Well 9: control. 
doi:1 0.1 371 /journal.pone.01 05404.g006 
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amount of (3-BuTx added to coat the plate wells in each round was 
fixed to 1 [Lg, and the amount of input ssDNA was gradually 
decreased (Table SI). BSA was used to block the wells in rounds 1- 
5 and then casein replaced BSA in rounds 6-10. The switch of 
blocking agent could dramatically decrease aptamer pool binding 
to the plate surface [36], so negative selections using uncoated 
wells were not included in our selection process. During the 
selection, the selection progress was monitored by determining the 
binding capacity of the ssDNA pool to (3-BuTx via nitrocellulose 
filter binding assay. Nitrocellulose filter membrane possesses long 
pores that are nucleic acid permeable but retain proteins as small 
as 2 kDa by hydrophobic adsorption [46]. The ssDNA-protein 
complex will retain on the filter when passing through the 
membrane. As shown in Figure 1, with increasing numbers of 
selection rounds, the affinity of the DNA for the target increased 
accordingly, and totally ten rounds of selection were performed. 

Sequencing and structural analysis 

The elution product of the 10 th round selection was amplified 
with unlabeled forward and reverse primers. After purification, the 
dsDNA was ligated into pMD19-T clone vector and transformed 
to E. coli cells. Forty-five positive colonies were randomly selected 
for sequencing, and 37 sequences were returned. Sequence 
comparisons revealed that only four individual binding species 
remained in the population, and the sequences were designated as 
fJB-1, (5B-19, PB-20 and PB-32, respectively (Table 1). One of 
these individuals, PB-1, dominated the selection and constituted 
more than 70% of the clones, while three other sequences together 
represented 27% of the clones. The length of variable region of 
aptamer PB-20 (51 nt) was significantly longer than the 40 nt 
randomized region of ssDNA library. Since longer or shorter 
byproducts usually appeared during the PCR reactions, maybe 
some longer products, but not long enough to be eliminated by gel 
purification, existed in the final aptamer pool. 

A further analysis with these sequences revealed a few conserved 
short sequences between PB-1 and PB-20 (Table 1). The 
secondary structures of the four sequences were evaluated at 
25°C, 150 mM NaCl, and 2 mM MgCl 2 using the Mfold software 
(Figure S2). All predicted secondary structures have stem-loop 
architectures. More notably, aptamers PB-1 and PB-20 share 
similar types of stems and hairpin loops, and the conserved 
sequences are involved in the formation of stems and loops (Figure 
S2). 

Binding characterization of selected aptamers 

To evaluate the binding affinities of these sequences, fluores- 
cence anisotropy measurements were conducted for all the selected 
aptamers. The complete sequences of aptamers were synthesized 
and labeled with FAM at the 5 '-terminal, and then titrated with 
different concentrations of P-BuTx. When the fluorescendy 
labeled aptamer bound to the target protein, the molecular weight 
of the complex increased which hindered the rotation diffusion 
rate of the labeled fluorescein and increased polarization, resulting 
in a significant anisotropy change of the aptamer. Figure 2 shows 
the change in anisotropy of aptamers PB-1 and PB-20 versus P- 
BuTx concentration. The K d s were analyzed by using a single-site 
ligand binding model and those of PB-1 and PB-20 were found to 
be 65.9±7.8 nM and 83.8±4.3 nM, respectively. The K rf s of PB- 
19 and PB-32 were much higher. The K ( i of PB-19 was 
530±58 nM, while that of PB-32 was 995± 138 nM (Figure S3). 
It indicated that the affinity of PB-19 and PB-32 was much lower 
than that of PB-1 and PB-20. 

In order to compare the binding sites of the selected aptamers, 
competition assays were performed. Two nM of fluorescendy 



labeled aptamer (PB-1, PB-20, PB-19 or PB-32) was firsdy 
incubated with 288 nM of P-BuTx, and then titrated by unlabeled 
aptamer (3B- 1 . If the unlabeled (3B- 1 competes for the same 
binding site as labeled aptamers, the displacement of the aptamer 
by PB-1 will result in a decrease in anisotropy [47]. As shown in 
Figure 3, with increase amount of PB-1, a decrease of fluorescene 
anisotropy was observed. It indicated that labeled aptamers were 
displaced by competitor PB-1, suggesting that the binding sites of 
the four aptamers are the same. Control experiments were done 
by unlabeled random starting ssDNA library competed with 
labeled PB-1, and no apparent anisotropy changes were experi- 
enced (Figure 3). 

Binding capacity of aptamer (3B-1 against related 
proteins 

The aptamer (3B- 1 with highest affinity against P-BuTx was 
chosen for assessment of its binding capacity against related 
proteins. We compared the ability of PB-1 to bind P-BuTx, BSA, 
casein and ot-BuTx. BSA and casein were the proteins that used 
for blocking during selection process, and ot-BuTx is a postsynaptic 
neurotoxin from B. multicinctus venom. The results from 
nitrocellulose filter binding assay (Figure 4A) showed that about 
45% of PB-1 bound to P-BuTx, and less than 4% bound to other 
three proteins when 1 0 nM of PB- 1 incubated with 1 00 nM of P- 
BuTx, BSA, oc-BuTx and 2.1 u,g/mL of casein, respectively. The 
affinity of PB-1 to the four proteins was also measured via 
fluorescene anisotropy measurements. As shown in Figure 4B, 
when 2 nM of PB-1 incubated with 50 nM of P-BuTx and 10 fold 
excess of other three proteins, the anisotropy change caused by P- 
BuTx was significantly higher than those by other proteins. The 
results confirmed that the binding between PB- 1 and P-BuTx was 
not only highly affinitive, but also highly selective. 

Aptamer PB-1 specifically recognize B. multicinctus 
venom 

Snake venoms contain a highly complex mixture of proteins, 
enzymes, and various other substances with toxic and lethal 
properties. Moreover, some snake venoms are reported to contain 
enzymes that have DNase activity which can hydrolyze DNA 
[48,49]. To determine whether the high concentration of other 
proteins presented in the snake venoms would have an influence 
on the aptamer binding or hydrolyze aptamer DNA, the reaction 
between PB-1 and five snake venoms, B. multicinctus, B. fasciatus, 
Deinagkistrodon acutus, Gloydius brevicaudus and N. Naja Atra, 
was investigated by EMSA. B. multicinctus and B. fasciatus are 
from the same genus of family Elapidae, D. acutus and G. 
brevicaudus are from two different genera of family Viperidae, and 
N. Naja Atra is from genus Naja of family Elapidae. B. 
multicinctus venom contains P-BuTx, while the four other snake 
venoms do not. As shown in Figure 5, aptamer/protein complex 
can be seen on the gel when 100 nM of (3B- 1 mixed with 240 nM 
and 480 nM of P-BuTx, or with 15 and 30 ug/mL of B. 
multicinctus venom. While the concentration of P-BuTx and B. 
multicinctus venom increased to 960 nM and 50 u,g/mL, respec- 
tively, the aptamer DNA stuck in the sample well completely. 
When N. Naja Atra venom was tested at the concentration of 
50 u,g/mL, aptamer DNA was totally hydrolyzed. However, the 
other three snake venoms did not obviously hydrolyze aptamer 
PB-1 or bind to it when tested at the same concentration, as there 
were only free DNA bands on the gel, and the lighteness of the 
bands were almost the same as the control PB-1 band (Figure 5). 
The results here showed that the four snake venoms, B. 
multicinctus, B. fasciatus, D. acutus and G. brevicaudus venom, 
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did not apparendy hydrolyze aptamer (3B- 1 , and (3B- 1 only bound 
to B. multicinctus venom but not other structurally or functionally 
related three snake venoms, even from the same genus. 

To further confirm the specificity of (3B- 1 against (3-BuTx and 
B. multicinctus venom, ELAA was performed using the above 
venoms excluding N. Naja Atra venom and proteins P-BuTx, 
BSA, casein and a-BuTx. Snake venoms and proteins of 10 [Lg/ 
mL were coated onto ELISA plate wells and incubated with 10 
pM biotin-labeled aptamer PB- 1 , respectively. Then HRP-labeled 
Streptavidin was added and substrate OPD was used to achieve 
colorimetric responses. Control well was coated with (5-BuTx but 
no aptamer was added. As shown in Figure 6, the color of wells 
coated with (5-BuTx and B. multicinctus venom turned to be 
brown, while the wells of other proteins and snake venoms, as well 
as the control well, underwent no obvious optical change. Without 
the aid of a detection device, one can apparendy distinguish P- 
BuTx from BSA, casein and a-BuTx, and B. multicinctus venom 
from three tested snake venoms. Furthermore, the color of (5-BuTx 
well was deeper than that of B. multicinctus venom well. This was 
because that only part of the B. multicinctus venom components 
was P-BuTx. The results proved that PB- 1 bound only to P-BuTx 
and P-BuTx-containing B. multicinctus venom, but not to other 
proteins or snake venoms. It also confirmed that the binding 
between (3B- 1 and P-BuTx was not interfered with by other 
ingredients presented in the four species of snake venoms. 

In conclusion, we have used plate SELEX to successfully isolate 
DNA aptamers for P-BuTx with nanomolar affinity in 10 rounds 
of selection. Based on the results of the present study, it 
demonstrates that the aptamer PB- 1 specifically binds to P-BuTx, 
and can discriminate B. multicinctus venom from different species 
of snake venoms. These results suggest that aptamer PB-1 could be 

References 

1. Warrcll DA (2010) Snake bite. Lancet 375: 77-88. 

2. Chippaux JP (1998) Snake-bites: appraisal of the global situation. Bull World 
Health Organ 76: 515-524. 

3. Kasturiratne A, Wickremasinghe AR, de Silva N, Gunawardena NK, 
Pathmeswaran A, et al. (2008) The Global Burden of Snakebite: A Literature 
Analysis and Modelling Based on Regional Estimates of Envenoming and 
Deaths. PLoS Med 5: e218. 

4. Gutierrez JM, Warrcll DA, Williams DJ, Jensen S, Brown N, ct al. (2013) The 
Need for Full Integration of Snakebite Envenoming within a Global Strategy to 
Combat the Neglected Tropical Diseases: The Way Forward. PLoS Ncgl Trop 
Dis 7: c2162. 

5. Lalloo DG, Theakston RDG (2003) Snake Antivcnoms: Antivenoms. Clin 
Toxicol 41: 277-290. 

6. Gutierrez JM, Theakston RDG, Warrell DA (2006) Confronting the Neglected 
Problem of Snake Bite Envenoming: The Need for a Global Partnership. PLoS 
Med 3: el 50. 

7. WHO (2010) WHO Guidelines for the production, control and regulation of 
snake antivenom immunoglobulins 2010 World Health Organization Geneva 

8. GaoJF, Wang J, Qu YF, Ma XM, Ji X (2013) Immunorcactivity between 
venoms and commercial antiserums in four Chinese snakes and venom 
identification by species-specific antibody. J Immunol Methods 387: 211-218. 

9. Isbistcr GK, O'Lcary MA, HaganJ, Nichols K, Jacoby T, ct al. (2010) Cross- 
neutralisation of Australian brown snake, taipan and death adder venoms by 
monovalent antibodies. Vaccine 28: 798-802. 

10. Dong LV, Quycn LK, Eng KH, Gopalakrishnakone P (2003) Immunogcnicity 
of venoms from four common snakes in the South of Vietnam and development 
of ELISA kit for venom detection. J Immunol Methods 282: 13-31. 

11. Selvanayagam ZE, Gopalakrishnakone P (1999) Tests for detection of snake 
venoms, toxins and venom antibodies: review on recent trends (1987-1997). 
Toxicon 37: 565-586. 

12. Keefe AD, Pai S, Ellington A (2010) Aptamers as therapeutics. Nat Rev Drug 
Discov 9: 5.37-550. 

13. Hong P, Li W, Li J (2012) Applications of aptasensors in clinical diagnostics. 
Sensors 12: 1181-1193. 

14. Ellington AD, SzostakJW (1990) In vitro selection of RNA that bind specific 
ligands. Nature 346: 818-822. 

15. Tucrk C, Gold L (1990) Systematic evolution of ligands by exponential 
enrichment: RNA ligands to bacteriophage T4 DNA polymerase. Science 249: 
505-510. 



useful in the development of diagnostic assays for snakebites. As 
aptamers are susceptible to nuclease-mediated cleavage, they are 
not stable in biofluid samples. Further chemical modifications of 
the selected aptamer should be made to improve its biostability. 

Supporting Information 

Figure SI Representative pilot PCR of the 10 th round of 
selection. Electrophoresis on 3% agarose after 18-26 cycles of 
PCR amplification. Lane 1 : 20 bp DNA Ladder (TaKaRa); Lanes 
2-5: 18, 20, 22, 24 and 26 cycles of PCR. The products of 20 
cycles of PCR amplification were relatively specific fragments. 
(TIF) 

Figure S2 The secondary structures of aptamers PB-1, 
PB-19, PB-20 and PB-32. 

(TIF) 

Figure S3 K,i measurements of PB-19 and pB-32 by 
fluorescence anisotropy. A) PB-19, K d = 530±58 nM; B) PB- 
32, K d = 995±138 nM. 
(TIF) 

Table SI Selection parameters of each SELEX round. 

(DOC) 

Author Contributions 

Conceived and designed the experiments: FPY QSF Jinglin Wang. 
Performed the experiments: FPY XW YH XLT Jie Wang. Analyzed the 
data: FPY YZ QSF TZ WWX. Contributed reagents/materials/analysis 
tools: FPY QSF Jinglin Wang. Contributed to the writing of the 
manuscript: FPY QSF Jinglin Wang. 



16. McKeague M, DcRosa MC (2012) Challenges and opportunities for small 
molecule aptamer development. J Nucleic Acids 2012. 

17. Bunka DH, Stockley PG (2006) Aptamers come of agc-at last. Nat Rev 
Microbiol 4: 588-596. 

18. Tan W, Donovan MJ, Jiang J (2013) Aptamers from cell-based selection for 
bioanalytical applications. Chem Rev 113: 2842-2862. 

19. Hung HT, Hojer J, Thi Du N (2009) Clinical features of 60 consecutive ICU- 
treated patients envenomed by Bungarus multicinctus. Southeast Asian J Trop 
Med Public Health 40: 518. 

20. Chan JC, Cockram CS, Buckley T, Young K, Kay R, et al. (1995) Evenoming 
by Bungarus multicinctus (many-banded krait) in Hong Kong. J Trop Med Hyg 
98: 457-460. 

21. Pc T, Myint T, Htut A, Htut T, Myint AA, ct al. ( 1 997) Envenoming by Chinese 
krait (Bungarus multicinctus) and banded krait (B. fasciatus) in Myanmar. 
Trans R Soc Trop Med Hyg 91: 686-688. 

22. Qin GP (1998) China poisonous snake research. Guangxi Science and 
Technology Press, Nanning. pp. 708-753 [in Chinese]. 

23. Chu NS (2005) Contribution of a snake venom toxin to myasthenia gravis: the 
discovery of ct-bungarotoxin in Taiwan. J Hist Neurosci 14: 138—148. 

24. Chang CC (1999) Looking back on the discovery of a-bungarotoxin. J Biomed 
Sci 6: 368-375. 

25. Kondo K, Toda H, Narita K, Lee CY (1982) Amino Acid Sequences of Three 
P-Bungarotoxins (P3-, [}4-, and (35-Bungarotoxins) from Bungarus multicinctus 
Venom. Amino Acid Substitutions in the A Chains. J Biochcm 91: 1531-1548. 

26. Kelly RB, Brown FR (1974) Biochemical and physiological properties of a 
purified snake venom neurotoxin which acts presynaptically. J Neurobiol 5: 
135-150. 

27. Chuang LY, Lin SR, Chang SF, Chang CC (1989) Preparation and 
characterization of monoclonal antibody specific for a-bungarotoxin and 
localization of the epitope. Toxicon 27: 211-219. 

28. Strong PN, Gocrkc J, Oberg SG, Kelly RB (1976) bcta-Bungarotoxin, a pre- 
synaptic toxin with enzymatic activity. Proc Natl Acad Sci 73: 178-182. 

29. Kwong PD, McDonald NQ, Siglcr PB Hcndrickson WA (1995) Structure of 
beta2-bungarotoxin: potassium channel binding by Kunitz modules and targeted 
phospholipase action. Structure 3: 1109-1119. 

30. Ueno E, Rosenberg P (1996) Mechanism of action of P-bungarotoxin, a 
presynaptically acting phospholipase A 2 neurotoxin: Its effect on protein 
phosphorylation in rat brain synaptosomes. Toxicon 34: 1219-1227. 



PLOS ONE | www.plosone.org 



7 



August 2014 | Volume 9 | Issue 8 | e105404 



Snake Venom Recognition by Aptamer 



31. Lauridscn LH, Vccdu RN (2012) Nucleic acid aptamcrs against biotoxins: a new 
paradigm toward the treatment and diagnostic approach. Nucleic Acid Ther 22: 
371-379. 

32. Tok JBH, Fischer NO (2008) Single microbead SELEX for efficient ssDNA 
aptamer generation against botulinum neurotoxin. Chcm Commun 16: 1883- 
1885. 

33. Cruz-Aguado J A, Penner G (2008) Determination of ochratoxin A with a DNA 
aptamer. J Agric Food Chcm 56: 10456-10461. 

34. Hcssclbcrth JR, Miller D, Robertas J, Ellington AD (2000) In vitro selection of 
RNA molecules that inhibit the activity of ricin A-chain. J Biol Ghem 275: 
4937-4942. 

35. Lauridscn LH, Shamaileh HA, Edwards SL, Taran E, Vcedu RN (2012) Rapid 
One-Step Selection Method for Generating Nucleic Acid Aptamcrs: Develop- 
ment of a DNA Aptamer against a-Bungarotoxin. PloS one 7: c41702. 

36. Hicke BJ, Marion C, Chang YF, Gould T, Lynott CK, et al. (2001) Tenascin-C 
aptamcrs arc generated using tumor cells and purified protein. J Biol Chcm 276: 
48644-48654. 

37. Lou X, Qian J, Xiao Y, Vicl L, Gcrdon AE, ct al. (2009) Micromagnetic 
selection of aptamcrs in microfluidic channels. Proc Natl Acad Sci 106: 2989- 
2994. 

38. Svobodova M, Pinto A, Nadal P, O'Sullivan CK (2012) Comparison of different 
methods for generation of single-stranded DNA for SELEX processes. Anal 
Bioanal Chem 404: 835-842. 

39. McEntee K, Weinstock GM, Lehman IR (1980) rccA protein-catalyzed strand 
assimilation: stimulation by Escherichia coli single-stranded DNA-binding 
protein. Proc Natl Acad Sci 77: 857-861. 



40. Morris KN.Jcnscn KB, Julin CM, Weil M, Gold L (1998) High affinity ligands 
from in vitro selection: complex targets. Proc Natl Acad Sci 95: 2902-2907. 

41. Levesquc D, Bcaudoin J, Roy S, Perreault J (2007) In vitro selection and 
characterization of RNA aptamcrs binding thyroxine hormone. Biochem J 403: 
129-138. 

42. Jocng CB, NiaziJH, Lee SJ, Gu MB (2009) ssDNA aptamcrs that recognize 
diclofenac and 2-anilinophcnylacetie acid. Bioorg Med Chcm 17: 5380-5387. 

43. Mendonsa SD, Bowser MT (2004) In vitro evolution of functional DNA using 
capillary electrophoresis. J Am Chem Soc 126: 20—2 1 . 

44. Hui Y, Shan L, Lin-fu Z, Jian-hua Z (2007) Selection of DNA aptamcrs against 
DC-SIGN protein. Mol Cell Biochem 306: 71-77. 

45. Hu P, Liu Z, Tian R, Ren H, Wang X, et al. (2013) Selection and identification 
of a DNA aptamer that mimics saxitoxin in antibody binding. J Agric Food 
Chem 61: 3533-3541. 

46. Ryan PC, Lu M, Draper DE (1991) Recognition of the highly conserved 
GTPase center of 23 S ribosomal RNA by ribosomal protein LI 1 and the 
antibiotic thiostrcpton. J Mol Biol 221: 1257-1268. 

47. Li W, Wang K, Tan W, Ma C, Yang X (2007) Aptamcr-based analysis of 
angiogenin by fluorescence anisotropy. Analyst 132: 107-113. 

48. Sittcnfeld A, Raventos H, Cruz R, Gutierrez JM (1991) DNase activity in Costa 
Rican crotaline snake venoms: quantification of activity and identification of 
cleetrophoretic variants. Toxicon, 29: 1213—1224. 

49. dc Roodt AR, Litwin S, Angel SO (2003) Hydrolysis of DNA by 17 snake 
venoms. Comparative Biochemistry and Physiology Part C: Toxicology & 
Pharmacology 135: 469-479. 



PLOS ONE | www.plosone.org 



8 



August 2014 | Volume 9 | Issue 8 | e105404 



